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ABSTRACT 

We have conducted an Ha survey of 38 0.77 < z < 1 galaxies over ~ lOOarcmin 2 
of the Hubble Deep Field North and Flanking Fields, to determine star formation rates 
(SFRs), with the near-infrared multi-object spectrograph CIRPASS on the WHT. This 
represents the first successful application of this technique to observing high redshift 
galaxies. Stacking the spectra in the rest-frame to infer a total SFR for the field, we 
find a lower limit (uncorrected for dust reddening) on the star formation rate density 
at redshift z = 1 of 0.04M Q yr _1 Mpc~ 3 . This implies rapid evolution in the star 
formation rate density from z — to z — 1 which is proportional to (1 + z) . 

Key words: instrumentation: spectrographs - galaxies: evolution - galaxies: forma- 
tion - galaxies: high-redshift - infrared: general 



1 INTRODUCTION 

The quest for a global star formation history of the universe, 
as a key element in understanding galaxy assembly, com- 
menced around 10 years ago. Madau et al. (1996) combined 
the results from the Canada-France Redshift Survey (CFRS) 
work (Lilly et al., 1995, 1996), Steidel's Lyman break galax- 
ies (LBGs) at z w 3 (Steidel et al. 1996a,b) and a population 
of faint galaxies with photometric redshifts 2 < z < 4.5 in 
the Hubble Deep Field (HDF), to give the first global view 
of the history of star formation in the Universe. Diagrams 
in a similar vein were produced around the same time by 
Lilly et al. (1996), at z < 1 and Pei & Fall (1995), based 
on damped Lyman alpha systems. Since those major works, 
there has been a volley of attempts to populate the now 
famous 'Madau-Lilly' diagram with data points at different 
epochs, using observational studies in a variety of wavebands 
to infer the (volume-averaged) star formation history. 

To date, many quantitative attempts to measure the 
global star formation history have been based on optical 
measurements and have thus suffered from having to use dif- 
ferent indicators of star formation in various redshift bins, 
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redshifted into the optical. These various indicators not only 
have uncertain relative calibration but are also affected dif- 
ferently by dust extinction. Commonly used star formation 
rate (SFR) indicators are ultra-violet (UV) continuum lu- 
minosity, which can be heavily dust extincted, and nebular 
emission lines such as Ha and [Oil] (the latter of which is 
strongly dependent on metallicity and ionisation state). A 
thorough review of these different techniques and their rel- 
ative calibrations is given by Kennicutt (1998). 

Longer wavelength estimators relying on Far-infrared 
(FIR) or radio luminosity are insensitive to dust obscura- 
tion yet have their own caveats. For example, radio lumi- 
nosity can also be produced by AGN, therefore this method 
requires cross-correlating with optical spectroscopic features 
to eliminate AGN from the sample. Sub-millimetre estimates 
rely heavily on assumptions about the dust temperature 
(e.g. Hughes et al. 1998) and a significant fraction of sub- 
mm galaxies also turn out to harbour AGN (Chapman et al., 
2005). 

The Universe at redshift z ~ 1 — 2 is believed to be one 
of the most active epochs in galaxy formation and evolu- 
tion. Indeed, it is inferred to be the epoch at which large 
elliptical and spiral galaxies are assembled and therefore 
may also be the period of peak star formation in the Uni- 
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verse. Yet the Universe at this epoch is still not well stud- 
ied nor understood. Observations have long been hampered 
by the difficulties of observing objects at these redshifts in 
the visible wavebands. At redshift z ~ 1, key diagnostic 
spectral features are redshifted out of the optical into the 
near-infrared, which is a difficult regime to work in, and the 
rest-UV Lyman— a line is not accessible in the optical un- 
til around redshifts z ~ 2.5. The redshift range z ~ 1 — 2 
has hence traditionally been dubbed the spectroscopic 'red- 
shift desert', although there has been some recent progress 
in exploring this epoch through extending colour-selection 
techniques to star-forming galaxies at z ~ 1.5 — 2 (e.g. Adel- 
berger et al. 2004 , Steidel et al. 2004). The Gemini Deep 
Deep Survey (Abraham et al., 2004) and the K20 survey 
(Cimatti et al., 2002) have also identified some passively 
evolving galaxies at these epochs. 

There is certainly evidence that the star formation rate 
was much higher in the recent past (z > 0.5), compared with 
the current epoch (e.g. Lilly et al. 1996, Tresse et al. 2002, 
Hippelein et al. 2003). However, it is still unclear whether 
at redshifts of one and beyond the star formation density 
plateaus, or declines or perhaps continues to increase. 

Even in the most up-to-date versions of the Madau-Lilly 
diagram (e.g. Hopkins 2004 ) there is still almost an order of 
magnitude discrepancy between different indicators, at z ~ 
1. In order to obtain a self-consistent global picture of the 
SF history of the Universe, the same robust indicator needs 
to be employed across all redshift bins. Researchers have 
begun to converge on the Ha emission line as such a robust 
indicator. It is a good tracer of the instantaneous SFR as it is 
directly proportional to the ionising UV Lyman continuum 
radiation from the most short-lived (tMS < 20 Myr) and 
massive (> IOMq) stars. Ha has been widely used in surveys 
at low redshift (e.g. Gallego et al. 1995). It is particularly 
suitable as it is relatively immune to metallicity effects and is 
much less susceptible to extinction by dust than the rest-UV 
continuum and Lyman- a (which is also selectively quenched 
through resonant scattering). The calibration required to 
derive SFRs from Ha assumes that no Lyman continuum 
photons escape the galaxy, which is almost true - at most 
the escape fraction may be around 5% (see discussion in 
Glazebrook et al. 1999 and references therein). 

Although the models also depend on star formation his- 
tory, galaxy age and metallicities, it is the IMF to which 
the SFR calibration is most sensitive as the Ha emission 
is produced almost entirely by the most massive stars. The 
Kennicutt (1998) conversion used in this paper is calculated 
for Case B recombination at T e = 10, 000K (e.g. Osterbrock 
1989) and assumes solar metallicities and a Salpeter (1955) 
IMF 1 with mass cut-offs 0.1 and 1OOM0. The conversion is 
given by: 



SFR(M Q yr" 1 ) = 7.9 x 10~ 42 L(Ha) (ergs" 1 ) 



(1) 



At z ~ 1 Ha is shifted into the near-infrared and un- 
til recently near-infrared spectroscopy has been restricted to 
long-slit work and building samples using single object spec- 
troscopy is inefficient in terms of telescope time. Further- 
more, the small statistical samples obtained result in large 



A Scalo (1986) IMF yields star formation rates a factor of three 
times higher for the same Ha flux. 



uncertainties in the global properties of galaxies at z ~ 1. 
For example, Glazebrook et al. (1999) targeted a sample of 
13 z ~ 1 field galaxies in the CFRS and found star forma- 
tion rates 2-3 times higher than those inferred from the UV 
continuum. Tresse et al. (2002) have a larger Ha sample of 
33 field galaxies at 0.5 < 2 < 1.1 and have built the Ha 
luminosity function. Comparing this to their lower redshift 
sample (Tresse & Maddox 1998) and other work in the liter- 
ature, they find a strong rise in star formation rate - a factor 
~ 12 from z = 0.2 to z = 1.3. Erb et al. (2003) obtained a 
small sample (~ 16 galaxies) of Ha measurements at z ~ 2 
for the purposes of dynamical rotation curves but did not 
derive the SFR density. 

There has been successful ' mult i- object' grism spec- 
troscopy through slitless surveys in the near-infrared from 
space, but such an approach has poor sensitivity because 
of the high near- infrared background due to zodiacal light. 
In the HST/NICMOS survey of McCarthy et al. (1999) 33 
emission line galaxies were identified between 0.7 < z < 1.9. 
From this data set Yan et al. (1999) derived the Ha lumi- 
nosity function and therefore the SFR density, again finding 
a result which is an order of magnitude higher than the lo- 
cal Universe. Hopkins et al. (2000) repeated this experiment 
with an independent data set of 37 Ha emitting galaxies in 
~ 4.4arcmin 2 , over a redshift range z — 0.7— 1.8 and found 
a consistent result. 

It is only now that true multi-object infrared spec- 
troscopy is possible from the ground, using either a slitmask 
approach (e.g. IRIS 2 on AAT and FLAMINGOS on Gemini 
and KPNO) or a fibre-fed spectrograph. 

We have used the Cambridge InfraRed PAnoramic Sur- 
vey Spectrograph (CIRPASS; Parry et al. 2000) in multi- 
object mode with the aim of addressing the true star for- 
mation history of the Universe at redshifts z = 0.7 — 1.5, 
through Ha measurements of a large sample of galaxies. 

CIRPASS can operate with 150 fibres with the ability to 
simultaneously observe up to 75 targets (in object/sky pairs) 
and demonstrates a powerful new technique for studying 
distant galaxies. Initial results from this project have been 
presented in Doherty et al. 2004. This is the first successful 
example of near-infrared multi-object spectroscopy of high 
redshift galaxies. Here we present results from this data set 
on the statistical measurement of the volume emissivity of 
Ha at redshift one, and hence the total star formation rate 
density of the universe at this epoch. 

The format of the paper is as follows. In Section 2 wc 
give a brief overview of the instrument and the observations. 
A full description of the instrument set-up and data reduc- 
tion is given in Doherty et al. 2004. We discuss the sample 
selection and completeness in some detail in Section 3. Sec- 
tion 4 covers the procedures behind stacking the spectra and 
calculating the total star formation rate of our sample and 
in Section 5 we discuss the necessary corrections to estimate 
a total volume-averaged SFR at redshift one, from the com- 
bined SFR of our sample. The reader familiar with this field 
and interested in our result in the context of evolution in 
the SFR density over cosmic time, may wish to skip directly 
to Section 6, where we compare our measurement to other 
Ha results from 2 = — 2. 

In this paper we adopt the standard "concordance" 
cosmology of Qj» = 0.3, £Ia — 0.7, and use hro = 
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#o/'<'Okms- 1 Mpcr 1 . AB magnitudes (Oke & Gunn, 1983) 
are used throughout. 



(see Section 4.2). We are in a Poisson regime as multiple 
readouts were taken to reduce the readnoise to a negligible 
level. 



2 OBSERVATIONS AND DATA REDUCTION 

CIRPASS is a near-IR fibre-fed spectrograph operating be- 
tween 0.9 and 1.67/im. The upper cut-off is set by a blocking 
filter which reduces the thermal background. CIRPASS can 
operate in one of two modes - with an Integral Field Unit 
or in multi-object mode (Parry et al., 2000). The CIRPASS 
Multi-Object Spectrograph (CIRPASS-MOS) was used at 
the Cassegrain focus of the 4.2m WHT in La Palma, to ob- 
serve 62 objects at a time in the Hubble Deep Field North 
(HDFN; Williams et al. (1996)). The fibre size corresponds 
to ~1.1" at the WHT (which is approximately comparable 
to the expected seeing convolved with typical galaxy pro- 
files, at least for compact galaxies). However, there are some 
galaxies in our sample which are extended and closer to ~2" 
across (as evident from HST imaging in the HDFN) and for 
these some light is therefore lost. At the WHT the fibres are 
deployable over an unvignetted field of 15' diameter, with 
a minimum separation of 12.6". A Hawaii 2K detector was 
used and a grating of 831 lines/mm, producing a dispersion 
of 0.95A pixel" 1 . The FWHM of each fibre extends over 2.7 
pixels in both the spatial and spectral domain. The wave- 
length coverage with the 2K detector was 1726A, covering 
most of the J— band, and the grating was tilted to set a cen- 
tral wavelength of A c « 1.2/im. We acquired the field using 
six guide bundles, each containing seven closely packed fi- 
bres, centred on bright stars. We thus positioned the plate 
to an accuracy better than 0.5", that is, less than half the 
fibre size. 

The observations were taken with the CIRPASS-MOS 
on the WHT, spanning a total of 6.5hours over two nights. 
These are summarised in Table 1. Precise details of the in- 
strument setup and data reduction process are given in Do- 
herty et al. (2004). 

Flux calibration was carried out to account for both 
the spectral response of the instrument and the atmospheric 
transparency. Observations of a standard star were used to 
calibrate the shape of the spectral response function (nor- 
malised to one in the middle of the detector) and the at- 
mospheric absorption. The flux normalisation was then ob- 
tained from a bright 2MASS star (J = 10.68) which was 
observed simultaneously with the targets i.e. with the same 
fibre plug plate. As the star and objects were observed si- 
multaneously, this method corrects for temporal changes in 
the seeing, which fluctuated considerably over the course 
of the observations. A second, fainter (J = 13.247) star was 
also observed but may have been vignetted as it was close to 
the edge of the plate and the relative fluxes are consistent 
to only 30%. Furthermore, the line fluxes of our detected 
objects vary between the two nights within this 30% error 
margin. This might suggest that there was a small rotational 
error in aligning the plate and we can thus only guarantee 
the flux calibration to within 30%. 

The gain of the 2K detector was measured to be 
4.5 e~ DN _1 where DN is data numbers, or counts. We kept 
track throughout the reduction process of the number elec- 
trons in the background noise, for Poisson statistics used in 
calculating the noise and hence the error on the final SFR 



3 SAMPLE SELECTION AND 
COMPLETENESS 

Our sample was selected from Cohen et al. (2000) - a very 
complete, magnitude limited, spectroscopic follow-up of the 
HDF-N. This catalogue also includes the 'Flanking Fields' 
(FF) - an approximately circular area of radius 8' centred 
on the HDF-N and imaged with WFPC2 for ~2 orbits in a 
single passband and the size of which is well matched to the 
CIRPASS-MOS field of view. 

The Cohen et al. (2000) catalogue is > 90% complete 
to R < 23 in the flanking fields and R < 24 in the HDF- 
N proper (AB magnitudes). In the catalogue there are 157 
galaxies in our range 0.7 < 2 < 1.0, but we selected only 
galaxies identified with emission lines (i.e. excluding absorp- 
tion line systems classified 'A'), giving 137. These 'A' sys- 
tems not only make up a small fraction (13%) of all targets in 
our redshift range but are likely to be either post-starburst 
or highly evolved galaxies with no current star formation. 
We therefore would not expect to detect Ha in these and 
chose not to allocate fibres to them, as they are unlikely to 
contribute much to the SFR density. 

We observed a random sub-sample of 62 of those 137 
galaxies, with perhaps only a small bias against close galaxy 
pairs, as the physical size of the fibres meant they had to 
be placed at least 12.6" apart on the sky. Where there was 
a position clash within that distance, the brighter source in 
R— magnitude was given priority. 

Sources which have redshifts placing Ha in the atmo- 
spheric absorption trough in the J— band at 1.10 — 1.16/im 
were excluded from the sample a posteriori, as strong atmo- 
spheric features create spurious signals in the stacked spec- 
trum (Section 4). This region is at the short-wavelength end 
of the array and taking this cut effectively increases our min- 
imum redshift to z = 0.768, decreasing the surveyed volume 
but increasing the confidence in the stacked detection. 

Our observed target list includes a small number of 
galaxies in magnitude bins where the catalogue complete- 
ness is low (i.e R > 23.5 in the FF and R > 24.0 in HDF). 
In order to derive a meaningful statistical result for the star 
formation in the field, we need a consistent sample to some 
limiting magnitude. We therefore take a cut at R — 23.5 
where there is ~84% completeness in the catalogue (with 
most of the incompleteness in the last 0.5 magnitude bin in 
the FFs). 

In summary, the sample adopted for the purposes of a 
statistical measurement of the SFR is restricted to 0.768 < 
z < 1.0 and 7? < 23.5 and numbers 38 objects (out of 91 pos- 
sible targets for observation). We therefore adjust the results 
by dividing by this incompleteness factor (38/91)=0.42. 

Guide stars were taken from the Guide Star Catalogue- 
II (Morrison et al. 2001). For consistency in the astrometry 
between our targets and alignment stars, we redetermined 
the coordinates of each target object using the GOODS vl.O 
images (Giavalisco & GOODS Team 2003). The typical off- 
set from the Cohen et al. (2000) coordinates was 0.5" (note 
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Table 1. Observations in the HDF-N. The data were taken over two nights with a total of 6.5 hours exposure time. The telescope was 
nodded by a distance 6.08 '(bcamswitch offset) between the object and sky fibres. 



Field 
name 


field 
centre 
(J2000) 


Date 
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no. 
objects 
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Figure 1. Histograms in R- magnitude (top) and Redshift (bot- 
tom) for all objects in the Cohen et al. (2000) catalogue in our 
redshift range (dashed line) and our observed sample (solid line). 
There seems no significant bias in magnitude or redshift space, 
particularly for our restricted sample with 0.768 < z < 1.0 and 
R < 23.5. 



that for a few extended galaxies our centre is also slightly 
different from the position listed in the GOODS catalogue). 

As a test for any bias in redshift or magnitude space we 
plot the histograms for our observed sub-sample versus all of 
the possible targets from the catalogue in our survey volume 
(Figure 1). There seems no significant bias in magnitude or 
redshift space, particularly for our restricted sample with 
0.768 < z < 1.0 and R < 23.5. 

Finally, there is a source of incompleteness introduced 
by taking a broad-band limit for the selection. Sources with 
large Ha equivalent widths (EqW(Ha)), i.e. sources which 
are luminous in Ha but too faint in the R— band to make 
the R < 23.5 magnitude cut (Figure 2), may be missed. 



Figure 2. He? v's R— magnitude diagram. The shaded regions 
show the high EqW(Ha) sources we would miss (in the Ha— R- 
mag plane) by taking a cut of R = 23.5 and the sources which 
we don't detect, but which contribute to the stacked flux (i.e. the 
statistical measure of the SFR, see Section 4). The sources with 
> 5<r Ha detections are plotted as triangles. 



We might reasonably expect there to be a population of 
sources with faint continuum magnitudes but large Ha flux, 
perhaps young, actively star forming galaxies with mini- 
mal evolved populations (which tend to dominate the longer 
wavelengths). Such systems would be missed in surveys like 
this which rely on magnitude- limited catalogues to pre-select 
targets with known spectroscopic redshifts. 

However, narrow-band imaging searches which select di- 
rectly on the emission line flux (in this case Ha) find few 
candidates and would suggest that the hypothetical popu- 
lation proposed above is extremely rare (if it exists at all) . 
For example, Willis & Courbin (2005) carried out a deep 
narrow-band search in the J— band, targeted at Ly— a at 
z ~ 9 to a limiting flux of 3 x 10 _18 erg s _1 cm -2 , over 
2.5 x 2.5arcmin 2 . The search is centred at 1.187/im and one 
might expect the by-product of Ha detections at z — 0.8. In 
fact they find only two sources with a narrow-band excess, 
one of which is likely Ha and the other H/3. 

Glazebrook et al. (2004) used a tunable filter with pre- 
cisely the aim of searching for line-emitters which would be 
missed in surveys such as ours. They observed three slices at 
7100, 8100 and 9100A (covering [Oil], H/3, Ha over the red- 
shift range 0.3-0.9), to a flux limit of 2 x 10~ 17 erg s _1 
cm~ 2 . They find no evidence for any new strongly line- 



emitting population of low luminosity objects and find an 
SFR density consistent with estimates from broad-band se- 
lected samples. 

For a typical detected galaxy (Rab ~ 22; Doherty et al. 
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2004) in our Ha survey, the limiting equivalent width at our 
5 a flux threshold of 1.9 x 10~ 19 erg cm~ 2 s is EW iest > 35 A, 
assuming a colour of (R — J)ab ~ 1 typical of star-forming 
galaxies at z « 1. For the broadband selection limit Rab — 
23.5, an EW closer to 140 A would be needed to detect Ha 
emission above our limiting flux. 



4 STAR FORMATION RATE FROM 
COMPOSITE SPECTRUM 

Confident that we have observed a random sub-sample of 
the whole sample (as discussed in Section 3), we now stack 
the spectra in rest- frame wavelength to obtain statistically a 
firm lower limit on the total star formation rate per unit vol- 
ume in our surveyed field. We subsequently discuss various 
corrections that can be made to estimate the true star for- 
mation rate density (e.g. incompleteness due to luminosity 
bias, aperture corrections and reddening corrections) and 
finally place this measurement in the context of evolution 
from 2 ~ 2.5 to z ~ 0. 



4.1 Stacked spectrum 

Any remaining skyline residuals due to imperfect back- 
ground subtraction that are close in wavelength to the Ha 
line could skew a statistical measurement of Ha (although 
note that none of the Cohen redshifts place Ha within 4A 
of a skyline). We therefore firstly mask out and interpolate 
over skyline residuals, where the sky is brighter than some 
threshold value, in the observed frame. 

Using the Cohen et al. (2000) redshifts, which are 
mostly determined from [Oil], each spectrum is shifted to 
the rest-frame, interpolated onto a common wavelength grid, 
which is over-sampled to 0.5A per pixel and summed. 

The stacked spectra (using a straight sum with no 
weighting) are shown in Figure 3. The Ha signal in the 
stacked spectrum from the Night I data is broader than 
that of Night II, in particular due to the presence of a fea- 
ture just blue- ward of the Ha line (around ~6555A). It is 
possible that it may be due to additional Ha, perhaps from 
sources with slightly incorrect optical redshifts. We visually 
inspected the individual spectra contributing to the stack 
and investigated stacking them cumulatively, to ascertain 
whether this feature represents a real signal. The sky sub- 
traction for this night was not as effective as Night II, and 
from inspection we deduce that this may be an accumulation 
of several sky residuals (in spite of all efforts to minimise the 
presence of skylines). Furthermore, since there is no evidence 
for a repetition of this signal in Night II data (where the sky 
subtracted out in a much cleaner fashion), we conclude that 
it is spurious. 

There is no evidence of [Nil] in the stacked spectrum. 
However, the [Nil] emission is known to have a systemic 
velocity offset from Ha, which may have a smearing effect 
in the stack. 

If we exclude the objects with robustly detected Ha 
(> 5a) from the stack, there is still a significant combined 
Ha contribution from the nominally undetected objects - 
around half of the flux in the total stacked spectrum (Fig- 
ure 4). 
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Figure 3. Stacked spectrum for each night's observations in the 
HDF-N, excluding sources where Ha lands in the atmospheric 
trough and restricting to R < 23.5. 38 sources in total. The central 
wavelength of Ha (6563A; from optical redshifts) is marked with 
a dashed line. The expected positions of [Nil] are marked with 
dotted lines, but they arc undetected. 



4.2 Calculating the total stacked flux and 
luminosity 

For each object we measure the Ha flux and hence line lumi- 
nosity and SFR, given its redshift. The sum of these SFRs 
represents the total star formation rate for our sample of 
galaxies. 

Figure 5 shows the curve of growth for the total Ha flux 
for each night's stacked spectrum, as a function of the wave- 
length extraction half- width (in the rest-frame) . It is evident 
that beyond a half-width of 4A, the summed flux from night 
II is rapidly converging, whilst that of night I continues to 
rise, due to the gradual inclusion of the noise spike corre- 
sponding to the spurious feature seen in the stacked spec- 
trum. We therefore determined the Ha line flux for each 
object by integrating over a wavelength region of 6563 ± 4A 
rest-frame (i.e. within ~ 200km s _1 of the [Oil] redshift). 
This extraction width optimises the real signal included in 
the final sum whilst excluding additional noise. The average 
baseline width between zero power points for the > 5a de- 
tections is 13A in the observed frame (see example spectra 
in Doherty et al. 2004, Figure 1). Our extraction width used 
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Figure 4. Stacks for each night excluding those sources detected 
at > 5<x. 29 sources total. 
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Figure 5. Total Ho flux in the stack (in counts), as a function 
of the extraction half-width taken. 



here is slightly larger and therefore encompasses the small 
velocity offsets between the [Oil] and Htv emission. 

To combine the data from the two nights, which should 
improve the S/N, the weighted mean of the two luminosity 
measurements for each object is taken according to: 



EIUV-? (2) 

where X; are independent estimates of the data point with 
associated errors en, and the variance is given by: 

ff2= E£TW (3) 

This method effectively reduces the error, providing a 
cleaner measurement of the total Ha flux. The noise in the 
integrated flux for each spectrum is calculated using Poisson 
statistics from the background noise, rather than measuring 
the standard deviation between pixels in the final spectra, 
as the noise is correlated after repeated extractions, inter- 
polations and rebinnings. 

The conversion to total star formation rate for the sam- 
ple is then given by Kennicutt's (1998) calibration (Equa- 
tion 1). 

The fluxes, luminosities and SFRs are shown in Table 2. 
The resulting combined SFR is a factor of ~ 2 higher than 
the sum of the individual > 5u detections. There is clearly 
a significant contribution to the total Ha flux from objects 
which are nominally undetected in the individual spectra 
(Figure 4). This therefore becomes a powerful technique for 
obtaining the total SFR in a surveyed volume, particularly 
when one can quantify the sample completeness. 



5 GLOBAL STAR FORMATION RATE 
DENSITY 

The SFR quoted above (131 ± 39M yr -1 ) is a firm lower 
limit on the true SFR for the total sample, as only a sub- 
set were spectroscopically targeted. In the following sub- 
sections we apply a series of corrections, leading to an es- 
timate which is closer to the true SFR for the field, and 
hence derive the global star formation rate density (SFRD) 
at z ~ 1. We then compare this value to other Ha mea- 
surements of the SFRD in the context of evolution from a 
redshift of around z ~ 2.5 to the present-day. 

The first, simple, correction is to account for the in- 
completeness due to crowding (i.e. the fact that we did not 
observe all galaxies in the field down to the limiting mag- 
nitude, but only a random subsample which could fit on 
the plug- plate without fibre clashes). This was discussed in 
detail in Section 3 and amounts to dividing the results by 
0.42, giving a corrected SFR of 312 ± 93M yr" 1 . The re- 
maining corrections rely on various assumptions, which shall 
be clearly stated and discussed as we proceed. The steps, 
briefly, are as follows: 

1) Incompleteness correction for galaxies below the 
survey flux limit. This is an effect dependent on the inher- 
ent luminosity function - that is, we do not probe the entire 
survey volume down to the same luminosity limit. In reality 
we have an approximately uniform flux limit (determined 
by the sensitivity of the instrument) which means that the 
luminosity limit drops off as l/d% where d,L is the luminos- 
ity distance. In other words, more distant galaxies must be 
intrinsically more luminous in Ha in order to be detected 
within our survey. To account for this incompleteness and 
estimate the true SFR density, it is necessary to make some 
assumptions about the Ha luminosity function in our red- 
shift range. 
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Table 2. Ha fluxes, luminosities and SFRs for each nights' data. Measurements are listed for different subsets of the sample: i) all 
objects contributing to the stack, i.e. 38 objects with 0.768 < z < 1.0 and R < 23.5, ii) > 5<r detections only (9 objects), iii) sub-sample 
excluding the > 5cr detections (29 objects). 



Night I 


Sub- set 


# 


Ho flux 
counts s~ x 


S/N 


Ha flux 

(xlO- 15 ) 

erg s _1 cm - 2 A — * 


(xlO 42 ) 
erg s _1 


SFR. 

M Q /yr 


i) 
iii) 


38 
9 
29 


1.56±0.07 
0.79±0.03 

0.77±0.07 


21.4 
24.8 
11.7 


5.07±0.26 
2.61±0.11 
2.46±0.23 


18.28±0.86 
8.94±0.36 
9.34±0.78 


144.4±6.8 
70.6±2.9 
73.8±6.2 


Night II 


i) 
ii) 
iii) 


38 
9 
29 


1.47 ± 0.08 
0.67 ± 0.04 
0.8 ± 0.07 


17.5 
16.3 
10.9 


4.29 ± 0.25 
1.90 ± 0.12 
2.39 ± 0.22 


14.41± 0.84 
6.46 ± 0.39 
7.95 ±0.75 


113.8± 6.7 
51.1 ± 3.1 
62.8 ± 5.9 



2) Aperture corrections. The fibre size on the sky is 
only 1.1" in diameter, which is around the size of the com- 
pact galaxies, convolved with the seeing (seen in HST imag- 
ing, e.g. Figure 2 in Doherty et al. 2004) . However, for more 
extended galaxies we may be missing a substantial fraction 
of the Ha flux, particularly if there are localised star-forming 
regions which fall outside the fibre aperture, which indeed 
appears to be the case, at least for some of the galaxies in 
the sample (see, for example, objects a), c) f) in Figure 2 of 
Doherty et al. 2004). There is no way of ascertaining aper- 
ture corrections for the Ha flux directly from our data set, 
instead we will attempt to infer aperture corrections based 
on B— band photometry (corresponding to the rest-frame 
UV). 

3) Reddening correction. As Ha is nominally unde- 
tected in most of our targets, it is not possible to de- 
rive actual reddening values for all sources. From the 
B— band luminosities measured for aperture corrections 
we determine the SFRr/y and hence using the average 
SFR(Ha)/SFR(UV) ratio, derive an average reddening cor- 
rection using two different extinction laws - a Milky Way 
type law and the Calzetti law. This is a useful exercise to 
place our SFRD measurement in context with other work 
- many authors have applied reddening corrections to their 
data before deriving SFRDs. 



5.1 Incompleteness and the Galaxy Luminosity 
Function 

The luminosity distribution of galaxies in the Universe is 
known to be well-characterised by a Schechter (1976) func- 
tion, in which the number of galaxies per unit volume, 
<j}(L)dL, in the luminosity bin L:L + dL, is given by: 



4>{L)dL = <j)*(L/L*) a exp(-L/L*)d(L/L*) 



(4) 



This function is parametrised by (j>* , L* and a, where <j>* 
is a number per unit volume, V defines the 'knee' of the 
function, and a gives the 'faint end slope', the slope of the 
luminosity function in the log(</>)— log(L) in the region L <C 
L* . Equation 4 integrates to an incomplete gamma function, 
which can be solved numerically to find the expected number 
density n of galaxies brighter than a certain luminosity limit 



n(L > Li im ) = / <\>{V)dL = fT{a + 1, x) (5) 

and the total luminosity in those galaxies: 



L{L > L Hm ) = 0*L*T(a + 2, x) 



(6) 



where x = L/L*. For a faint end slope —2 < a < — 1, 
the function diverges for the number counts but not the 
luminosity (i.e. there can be an infinite number galaxies but 
not infinite luminosity - many very low luminosity galaxies 
contribute almost nothing to the total luminosity). However, 
for a faint end slope a < —2 there would also be infinite 
luminosity. 

Knowing the luminosity function to some degree of ac- 
curacy allows us to correct for the incompleteness in surveys 
introduced by the luminosity-distance dependence. In prac- 
tice, with only 9 solid detections, we do not have enough 
data points to construct a credible Ha luminosity function. 
However, we can attempt to model the incompleteness by 
investigating a range of Ha luminosity functions based on 
previous work 2 , comparing our observed result with the pre- 
dicted number counts and SFRs of star-forming galaxies 
which we should expect to see above the flux limit, over 
our redshift range. 

Given various input luminosity functions, parametrised 
by L* , 0*and a, we loop over redshift and over luminosity 
and for each redshift slice compute the expected numbers 
and SFRs in each luminosity bin (using equations 5 and 6), 
thereby producing a predicted number density as a func- 
tion of limiting line flux. We compare to our 5a flux limit 
in a single average fibre, /h„ = 1.9 x 10 _16 ergs _1 cm~ 2 
and to the effective 5<r flux limit obtained by stacking 
up the 38 observed galaxies in our survey field, fu a — 
3.3 x 10 _17 ergs _1 cm~ 2 . Table 3 gives a summary of the 
result of this investigation. It lists the input luminosity func- 
tion parameters, the corresponding predicted number counts 
and SFRs we should expect to observe in our surveyed vol- 
ume at the above flux limits, and the total SFR (obtained 
by integrating down the luminosity function). This method 



2 converted to our ACDM cosmology. See Appendix for the 
method behind this. 
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allows us to infer what order of correction is needed to esti- 
mate the total SFR for our field. 

Yan et al. (1999) derived an Ha luminosity function for 
galaxies at 0.7 < z < 1.9 using slitless spectroscopy with the 
near- infrared camera and multi-object spectrometer (NIC- 
MOS) on the HST. In that survey they detect 33 emission 
line galaxies over 64arcmin 2 (McCarthy et al., 1999). How- 
ever, we note that this number seems very low if their lumi- 
nosity function is correct (see Table 3 where the predicted 
number counts are high). Their average 3a line flux limit is 
4.1 x 10 _lr ergs _1 cm -2 (McCarthy et al., 1999). Using their 
derived luminosity function (LF) to model the galaxy counts 
and SFR we would expect to see in our survey volume, we 
find that either the L* or <j>* is significantly too large to 
account for our observed numbers - it predicts around 3-4 
times the number that we actually detect (Table 3). The 
large value of L* implies that the majority of the star for- 
mation is in galaxies brighter than our 5a limit, meaning 
if their luminosity function is accurate and appropriate for 
z ~ 1 galaxies, we should have detected most of the objects 
we observed in Ha. Hopkins et al. (2000) carried out a simi- 
lar experiment with NICMOS grism spectroscopy and find a 
similar star formation rate density, albeit a very differently 
shaped luminosity function. It is difficult to reconcile our 
result with those works if we take their measurement to be 
representative of the redshift z ~ 1 Universe. 

Turning now to compare our results to the work of 
Tresse et al. (2002) 3 , who derive an LF from ~ 30 emis- 
sion line galaxies at 0.5 < z < 1.1 in the CFRS, we find a 
somewhat closer match, but their LF still overpredicts our 
number counts and SFR by ~ 2.5 (Table 3). While the L* 
seems to reproduce a closer match to the proportion of the 
SFR we find in our detected galaxies (compared to L* found 
by Yan et al. (1999)) their value of <f>* seems too large to fit 
our data. In fact, this discrepancy is accounted for by aper- 
ture losses. Our fibre size is 1.1" diameter but Tresse et al. 
2002 used a 2" slit, which is better matched to the average 
galaxy size. Once we apply aperture corrections (detailed in 
the following Section 5.2) our result is comparable. 



5.2 Aperture Corrections 

There is no fool-proof method to calculate aperture correc- 
tions for our Ha fluxes. We attempt this using the GOODS 
B— band imaging, which at redshift z ~ 1 samples the rest- 
frame UV. We perform photometry on these images using 
phot in IRAF, with the same centres as our fibre positions, 
and 1" apertures. The aperture correction is then the dif- 
ference between the resulting 1" aperture magnitudes and 
the total magnitudes listed in the GOODS vl.O catalogue 4 . 
We then apply the same fractional correction to the Ha flux. 
This relies on the assumption that the distribution of Ha di- 
rectly traces the UV, that there is no gradient in reddening 
across the galaxy. 



3 Tresse et al. (2002) have applied reddening corrections to their 
data before fitting a luminosity function. However, they cite an 
average correction of A v = 1 mag (a factor of 2.02 in H« lumi- 
nosity) and we use this value to artificially 'redden' their LF for 
the sake of comparison. 

4 Publicly available at http://www.stsci.edu/science/goods/ 
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Figure 6. B— band 1" aperture magnitudes versus the fractional 
correction necessary to multiply by to obtain total B— magnitudes 
for our observed sample. The galaxies detected in Ha are circled. 
Note that we detect all but one galaxy brighter than B = 24. 
The overlaid object identification numbers allow a comparison 
with Table 3.6 and Figure 3.9 in the PhD thesis of Dohcrty (Uni- 
versity of Cambridge, 2005), which display the B— band (rest-UV) 
properties of these objects. 



Figure 6 shows the B— band 1" aperture magnitudes 
versus the fractional correction necessary to obtain total 
B— magnitudes for our observed sample, with the galaxies 
detected in Ha marked with circles. We detect all but one 
galaxy brighter than B = 24, indicating that perhaps the 
reason we don't detect more galaxies in Ha is simply the 
limit of the instrument sensitivity. 

As a consistency check we also looked at the 1" aperture 
magnitudes from the GOODS catalogue. However, we have 
not used these in calculating aperture corrections as our cen- 
tres are very slightly offset from the GOODS centres, in the 
most extreme case by ~ 0.5" for very large/irregular galax- 
ies. This therefore can significantly affect the 1" magnitudes 
- we need to compare to the same positioning as our fibres 
in order to extrapolate the Ha flux correction. It does not, 
however, affect the total nragnitudes and we have therefore 
taken these directly from the GOODS vl.O catalogue. 

Applying the measured fractional aperture corrections 
to the Ha flux measures (or upper limits) leads to an over- 
all aperture correction of a factor of 2.4 for the stack. 
Adjusting the SFR and SFRD estimates by this amount 
then gives 746 ± 224M yr" 1 and an SFRD of 0.034 ± 
0.010 Mq yr _1 Mpc~ 3 (for a co-moving volume of 2.2 x 
10 4 Mpc 3 ), which is more in line with the measurement of 
Tresse et al. (2002). From Table 3 it can be deduced that 
the completeness correction for luminosity bias is most likely 
of the order of 15-20%, giving ~ 0.04 M yr -1 Mpc" 3 . 



5.3 Reddening corrections 

Ha emission is more robust to extinction than the UV con- 
tinuum but not completely immune. In the following sec- 
tion we firstly compare uncorrected values of the Ha derived 
SFRD. However, this assumes that the amount of dust ex- 
tinction does not evolve over cosmic time and also that it 
is independent of the SFR (i.e. involves the same scaling 
over all redshift bins) which is known to be untrue (star 
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Table 3. Top two rows show our observed numbers and SFRs to both the limiting flux in an average fibre and the effective limiting flux 
(per object) when stacking the spectra. The numbers given are for the sample of 38 objects, corrected by a factor 0.42 for completeness. 
The second row includes aperture corrections. The lower limit on the SFRD is found by dividing the SFR by the total co-moving volume 
surveyed (V cm = 2.2 X 10 4 Mpc 3 ). We compare our result with various luminosity functions from the literature, determining predicted 
number counts and SFRs to the same limiting flux, in our surveyed volume. References are: Y99: Yan et al. (1999), T02: Tresse et al. 
(2002), G95:Gallego et al. (1995), H00: Hopkins ct al. (2000; two sets of parameters arc given corresponding to upper and lower limits 
on the Ho luminosity function). Also tabulated are a few test functions used to investigate the effect of varying L* and <fi*. Although 
not a rigorous fit this is a useful comparison and indicates a correction of the order of 15 — 20% for incompleteness due the luminosity 
bias. (L* is in units of erg s — 1 , <f>* in Mpc -3 .) 



Ref. Luminosity fn params a 

log(L*) <fr* c 

(xl0~ 3 ) 



< z > limiting flux 

/(Ha) > 1.9 x 10- 16 

erg s _1 cm~ 2 

# /50.25' 2 SFR (>lim) 



limiting flux 
/(Ha) > 3.3 x 10- 17 
cm * 
SFR (>lim) 



total SFR 

M /yr 



erg s 
# /50.25' 2 



SFR density 
M Q /yr/Mpc 3 



This paper 




observed 




0.82 


21 


167 


90 


312 


- 


>0.014 


This paper 


with 


ap. corrections 


0.82 


21 


283 


90 


746 


- 


>0.034 


Y99 


42.82 


1.5 


-1.35 


1.3 


83 


1782 


257 


2225 


2428 


0.109 


T02 


41.97 


4.3 


-1.31 


0.74 


26 


261 


205 


680 


922 


0.041 


G95 


41.56 


1.6 


-1.3 


0.0225 


1 


9 


33 


72 


133 


0.006 


H00 


42.88 


0.77 


-1.6 


1.3 


66 


1274 


285 


1789 


2284 


0.102 


H00 


43.35 


0.088 


-1.86 


1.3 


35 


778 


193 


1125 


2324 


0.104 


Test 


41.97 


2.0 


-1.35 


- 


12 


121 


100 


325 


456 


0.02 


Test 


41.97 


1.8 


-1.35 


- 


11 


109 


90 


293 


410 


0.018 


Test 


42.00 


4.0 


-1.35 


- 


27 


279 


212 


709 


977 


0.044 


Test 


42.00 


.3.8 


-1.35 


- 


26 


265 


201 


674 


928 


0.042 



n 42.S5 , 



0* = 1.7xl0- 3 Mpc- 3 ,a = -1.35 (for H =50 km s" 1 Mpc" 1 , q 
(Hq =50 km s — 1 Mpc — *, go = 0.5) reddening corrected; G95:a = 
Mpc" 1 ) reddening corrected; H00: 0*=3.3 X 10~ 4 ,L*=10 4302 , a 
Mpc" 1 , q = 0.5). 



-,42.:S7 
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-1 

_ w ^ -3 

4 Mpc" 3 (H = 50km s" 1 
* = -1.6 (H = 75km s" 1 



4.07xlO~ 3 Mpc 



formation creates dust therefore extinction is quite heavily 
dependent on the amount of ongoing star formation). The 
final correction needed therefore, is for dust reddening. 

The usual method for correcting Hq is to compare the 
observed Balmer decrement (Ha/H/3 ratio) with its the- 
oretical value (=2.86; Osterbrock 1989). We did not ob- 
serve the H/3 line as our wavelength coverage was insuf- 
ficient, we instead work with the UV-continuum from the 
B— band magnitudes measured above, and use the observed 
SFR(Hq):SFR(UV) ratio to derive an estimate of the ex- 
tinction using both the Calzetti (1997) law and a Milky Way 
law, for comparison with other authors (Section 6). 

We calculated rest-frame UV (2400 A) luminosity den- 
sities and corresponding star formation rate using the con- 
version given in Kennicutt (1998): 




SFR(M yr" 



1.4 x 10~ 28 L„ 



(erg s Hz 



(7) 



which assumes continuous star formation over ~ 10 years, 
a fiat continuum in the UV (in /„) and the same IMF as 
used when deriving SFRs from the Hq flux (Equation 1). 

Figure 7 shows the SFRs for each galaxy calculated us- 
ing the UV luminosity density and the Hq flux: those cal- 
culated from the UV luminosity densities are a factor of 
1.98 lower on average than those from Hq. This is proba- 
bly due to the differential effect of dust extinction in the 
redshift one galaxies between A rcst ~ 2400A and 6563A. 
This is consistent with results obtained by Glazebrook et 
al. (1999), Tresse et al. (2002), and Yan et al. (1999) who 
all find SFR(Hq)/SFR(UV) ratios of around 2-3. 

There is no uniform correction that can be applied to 



Figure 7. Comparison of SFRs obtained from UV continuum 
flux at 2400A versus Ha flux for the individual galaxies. The SFRs 
derived from UV luminosity arc consistently underestimated. The 
filled circles are the robust detections, the open circles are greater 
than 3<r and the squares represent 3<r upper limits. The solid line 
has a gradient of 1.98 and represents the line of best fit to the 
data (using a least squares fit, through zero). The errors in the 
SFR(Ha) values reflect the 30% error on the flux calibration. We 
assume ~ 5% error in the UV flux, which is dominated by errors 
in positioning the fibres (i.e. mirrored in the exact positioning of 
the aperture for the B^g magnitude). The dotted line is the line 
of zero extinction i.e. where SFR(Ha)=SFR(UV). 
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Table 4. E(B — V), A v , An a and correction factor to convert 
SFR(Ha) to a total SFR, for the Milky Way and Calzetti laws 



E(B - V) A v 



correction factor 



Milky Way 
Calzetti 



0.156 
0.16 



0.5 0.41 
0.78 0.65 



1.46 
1.85 



the SFR estimated from UV flux in order to derive the true 
star formation rate: the amount of extinction varies in each 
object due to inherent differences in physical properties of 
the galaxies (e.g. Sullivan et al. 2001). The errors on our 
Hq fluxes are not great enough to account for the scatter 
in Figure 7, which instead is attributable to different dust 
extinctions in our galaxy sample. However, we can estimate 
the true star formation rate for the sample, using the average 
SFR(Hq)/SFR(UV) ratio if we assume that the Kennicutt 
(1998) relations for Ha and UV SFRs should yield the same 
intrinsic SFR in the absence of extinction. 

Table 4 gives the values of the extinctions in magnitudes 
A v , Ana, the colour excess E(B — V) and the conversion 
factor to convert SFR(Hq) to a total SFR, for the Milky 
Way (Seaton 1979) and Calzetti laws (Calzetti 1997), which 
yield our observed ratio SFR(Ha)/SFR(UV)=1.98. For the 
purposes of comparison in the Madau-Lilly diagram in the 
following section, we show reddening corrections using both 
extinction laws. 



6 MADAU-LILLY DIAGRAM: EVOLUTION OF 
THE STAR FORMATION RATE DENSITY 

Having finally derived a global SFRD for the redshift range 
0.768 < 2 < 1.0, corrected for incompleteness, fibre aperture 
losses and with an average reddening correction, we turn 
now to placing this measurement in the context of other 
work in the field, to address the question of evolution of the 
SFR over cosmic time. As discussed in the Introduction, the 
past decade has seen many attempts to correct and calibrate 
various SFR indicators, in an attempt to put measurements 
at different cosmic epochs on the same scale (the 'Madau- 
Lilly diagram') and by so doing trace the time evolution of 
the global star formation rate density. We focus on compar- 
ing our result only to other Ha measurements in the liter- 
ature, in order to take one consistent indicator across all 
redshift bins. With current infra-red technology Ha can be 
traced out to redshift z ~ 2.5 (AT— band) - i.e. ~ 80% of the 
age of the universe. The compilation of Ha SFR measure- 
ments from the literature is given in Table 5, converted to 
our cosmology (H =70kms~ 1 Mpc~ , Qm = 0.3, £7a = 0.7). 
These are then used to plot the evolution of the SFRD (un- 
corrected for reddening) in Figure 8. Some authors have cor- 
rected their data for reddening and for the purposes of com- 
parison with our result, following Tresse et al. (2002), we 
artificially 'redden' their result using the canonical A v — 1, 
(a factor 2.02 at Ha). We overplot two values from this work 



5 Note that the correction from A v to An a is small, but non- 
ncgligiblc. Ajj a = 0.82j4„ for a Milky Way law and 0.85Ai, for a 
Calzetti law. 



- the first, lower point is our measured SFRD with no cor- 
rections (other than for the completeness of the observed 
sample) which represents a firm lower limit on the SFRD 
in the redshift range 0.768-1.0 . The upper value includes 
corrections for incompleteness due to luminosity biasing and 
aperture corrections. The error bars reflect the uncertainty 
in the flux calibration, which is ~ 30%. 

There may be a small amount of contamination from 
Active Galactic Nuclei (AGN). We have not corrected for 
this, both as the effect is very minimal (well within the 30% 
systematic errors) and also as the other SFR estimates on 
the Madau-Lilly diagram in Figure 8 have not been cor- 
rected. Pascual et al. (2001) estimate that at low redshift 
AGN contribute about 10% by number and ~ 15% of the 
Ha luminosity density. Matching our target list of galaxies 
observed with C1RPASS to the Chandra 2Ms X-ray cata- 
logue (Alexander et al. 2003), we find six sources out of our 
62 observed which are likely AGN. Only three of these are 
included in the stacked spectrum from which we estimate 
the SFRD, (the others were excluded due to Ha falling in 
the water absorption region of the spectrum) the combined 
SFR of these three sources as inferred from their Ha flux is 
12.5Mq yr~\ or 9.5% of the total for our observed sample. 
This constitutes an upper limit to the AGN contamination 
as not all of the Ha flux is produced by the AGN, it is likely 
a mixture from both the AGN and star formation. 

As some authors have corrected their data for redden- 
ing, we make a separate comparison to the reddening cor- 
rected SFRDs in the literature in Figure 9. We use the aver- 
age reddening on our data as derived above - a factor 1.46 
for a Milky Way law or 1.85 for a Calzetti law, to compare 
our result. As the authors who derived the three highest 
redshift data points (Yan et al. 1999, Hopkins et al. 2000 
and Moorwood et al. 2000) have not applied reddening cor- 
rections we adjust their data points by A v — 1 (or a factor 
2.02, following Tresse et al. 2002). 

Our SFRD at 0.768 < z < 1.0 is consistent with the 
value derived by Tresse et al. (2002) over a similar redshift 
range, although our data do not support quite the same 
strength of reddening as they invoke. Our result is still a 
factor of ~2.5 lower than the results of Yan et al. (1999) 
and Hopkins et al. (2000). However, these works span a much 
greater redshift range, probing out toz~ 1.8. 

The evolution in luminosity density, or star formation 
rate density, of the Universe is often fit by a power law 
of the form (1 + z) n . Fitting a least-squares line in the 
iog(SFRD)-log(l+z) plane yields (1 + z) 3 ' 1 (Figure 8), 
which is strong evolution from z = 1 to the z — universe, 
and comparable to previous work. Tresse et al. (2002) found 
evolution in the star formation rate density as n = 4.1 from 
Ha (for an Ho = 50km s _1 Mpc -1 , qo = 0.5 cosmology), 
which was in agreement with the evolution in UV (2800A) 
luminosity density in the CFRS (Lilly et al. 1996). In fact, 
the change in cosmology softens the power law. Without re- 
fitting their data, using Gallego's value as a zero point and 
comparing to that of Tresse, we find that the cosmology con- 
version would change the power law ton ~ 3.5. We there- 
fore find a power law slightly shallower than those previous 
works, but our data nevertheless supports a steep decline in 
the star formation rate from z — 1 to z — 0. 
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Table 5. Values in the literature for Ha luminosity densities (L(Ha) in erg s — 1 Mpc -3 ) or Ha derived star formation rate densities. 
Values of L(Ha) have been converted to our cosmology, Hq = 70 km s — 1 Mpc" 1 , A = 0.7, Q = 0.3, using the mean redshift < z >. 
Where the authors have applied a reddening correction to their data, we use a value of A„=l mag to 'de-redden' their result for the no 
reddening column. Conversely, if no reddening correction has been applied in the original work, we redden the result by A„=l mag to 
obtain the values listed under the reddening column. Star formation rate densities have been derived using the converted L(Ha) value and 
the K98 conversion. References are G95:Gallego et al. (1995) , PG03:Perez-Gonzalez et al. (2003) , S00:Sullivan et al. (2000) ,T98:Tresse 
& Maddox (1998) , P01:Pascual et al. (2001) , T02:Tresse et al. (2002) , D06: this paper, H00:Hopkins et al. (2000) , Y99:Yan et al. 
(1999) , M00:Moorwood et al. (2000) . 
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Figure 8. SFRDs determined from Ha measurements only, with 
no reddening corrections. Circles are points taken from the liter- 
ature, converted to our cosmology (see Table 5). Overlaid is our 
lower limit (square) to the SFRD, and our estimate including lu- 
minosity bias and aperture corrections (triangle) . The dashed line 
represents a (1 + z) 3,1 power law. Our corrected point is consistent 
with previous determinations of evolution in the SFRD according 
to (z + l) 4 (in an Ho = 50 km s — 1 Mpc -1 , go = 0.5 cosmology) 
from z = — 1 (e.g. Tresse et al. 2002 and refs therein). 



Figure 9. As for Figure 8 but with reddening corrections. We 
have used the authors own reddening corrections where applica- 
ble. For the three highest redshift points, the authors have not 
given reddening corrections and we use a correction A v = 1.0, 
or a factor 2.02 at Ha, following Tresse et al. (2002). Our un- 
corrected data points are shown with a square and triangle, as 
above. Reddening corrected points are shown for the Milky Way 
law (filled star) and a Calzetti law (open star). 



7 SUMMARY 

We have carried out an Ha survey of redshift one galaxies 
in the HDF-N, using multi-object near infrared fibre spec- 
troscopy. We found 9 robust detections (> 5a), of which the 
brightest are presented in Doherty et al. (2004). The real 



power of this technique lies in stacking the signal to obtain 
a star formation rate estimate for the sample as a whole. 

We have derived a lower limit to the SFR of 312 
Moyr -1 (corrected for incompleteness due to observing a 
random sub-sample of all possible targets), or an SFRD 
of 0.014 M0 yr~ Mpc - . Correcting this for volume in- 
completeness due to the luminosity bias effect and aperture 
losses gives an SFRD of 0.04 M Q yr" 1 Mpc -3 . This is con- 
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sistent with the work of Tresse et al. (2002) over a similar 
redshift range. Yan et al. (1999) and Hopkins et al. (2000) 
find a result ~ 2.5 times higher, from space-based grism 
spectroscopy, but their work covers a much wider redshift 
range (z ~ 0.7— 1.9), with a mean < z >=1.3 and we there- 
fore do not expect it to be representative of z = 0.8 if the 
Ha luminosity density continues to rise at z > 1. 

Our SFRD estimate is consistent with steep evolution 
from z = to z = 1, described by a power law (1 + z) . 

We finally derived an average reddening correction of 
A„=0.5-0.78 magnitudes, or a factor 1.46-1.85 at Ha (de- 
pending on the extinction law), and place this in context 
with reddening corrected values from the literature (Fig- 
ure 9). However, this may be an underestimate of the red- 
dening as it was derived using a comparison to the UV lumi- 
nosities, and the UV luminosity can continue to evolve even 
after star formation has ceased leading to an additional un- 
certainty in the SFRuv- A caveat to this work which should 
be noted is that the result is obtained for one small field of 
view and thus is cosmic variance dependent. 

We have demonstrated the first successful application 
of multi-object fibre spectroscopy to observe high redshift 
galaxies. In spite of the many difficulties involved in search- 
ing for such weak signals, this approach constitutes a pow- 
erful technique and the success of our multi-object spec- 
troscopic survey bodes well for larger surveys with future 
instruments such as FMOS on Subaru (Lewis et al., 2003; 
Kimura et al., 2003) and EMIR on GTC (Garzon et al. 
2004). 
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APPENDIX 

The conversion of an Ha luminosity function from the orig- 
inal cosmology to that assumed throughout this paper is 
as follows (following Hopkins 2004). For a flat universe, the 
co-moving volume in a certain redshift rang e z\ < z < Z2 
is proportional to the difference of the co-moving distances 
(Dc) cubed, i.e. V cm ex D c (z2) A — D c (zi) 3 and luminosity is 
proportional to the co- moving distance squared, L oc D^. 
To convert a luminosity function we convert the luminosity 
and number density (i.e. volume) parameters separately, by 
evaluating those expressions at the appropriate redshift, for 
each cosmology. The ratio gives the conversion factor, i.e. 
L* is multiplied by the ratio of the squares of the co-moving 
distances and <fi* is divided by the ratio of the co-moving 
volume shells (between z\ and z-i). So to convert a luminos- 
ity function parametrised as L* , <f>* , a in a cosmology where 
the co-moving distance is D c , to L*' , cj>*' , a' with D' c , wc 
have 



L*' = V x Df/Dl 



and 



D c {z 2 f - D c (zi) 3 



(8) 



(9) 



r D c (z 2 )' 3 - D c (zi)' 3 

This conversion relies on the assumption that all galax- 
ies lie at the central redshift, a more rigourous method would 
be to correct the original data points and refit the Schechter 
function. However, the likely error associated with employ- 
ing this approximation method is ~ 10% (Hopkins 2004). 



